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A new family of simple analytic solutions of hydrodynam- 
ics is found for slowly expanding, rotationally symmetric fire- 
balls assuming an ideal gas equation of state. The tempera- 
ture profile is position independent only in the collisionless gas 
limit. The Zimanyi-Bondorf-Garpman solution and the Buda- 
Lund parameterization of expanding hydrodynamic particle 
sources are recovered as special cases. The results are applied 
to predict new features of proton correlations and spectra at 
1.93 AGeV Ni + Ni collisions. 

Introduction — Hydrodynamical models are applied 
in a wide range of physical problems from the evolu- 
tion of galactic systems through the explosion of super- 
novas to the description of particle spectra and corre- 
lations in heavy ion collisions in a broad energy range. 
Due to the nonlinearity of the hydrodynamical equations 
it is rather difficult to solve them in an analytic man- 
ner. In particle and heavy ion physics some observables 
like the single-particle spectra and the two-particle Bose- 
Einstein correlation functions are sensitive to the means 
and variances of the hydrodynamical flow and density 
profiles 0^]. This result implies that a hydrodynamical 
solution involving only the means and the variances of 
the density profile, flow profile and temperature profile 
is of great interest. Such a solution of the equations of 
non-relativistic (NR) hydrodynamics is presented here- 
with. Physical realization of the solution is possible both 
by heavy ion collisions in the 30 MeV A - 2 GeV A en- 
ergy domain ||-^| and in a more limited sense by the 
NR transversal flows that seem to develop in high energy 
hadron-hadron and nucleus- nucleus collisions |l],|6||j]]. 

The new family of exact solutions — Let us consider 
the NR hydrodynamical problem, as specified by the con- 
tinuity, Euler and energy equations: 



d t n + V(vn) = 0, 
9 t v + (vV)v = -(Vp)/(mn), 
d t e + V(ev) = -pVv, 



(1) 
(2) 
(3) 



where n denotes the number density of particles, v stands 
for the NR flow field, e for the NR energy density, p for 
the pressure and in the following the temperature field 
is denoted by T. We assume a NR ideal gas equation of 
state, 



P = nT, 

3 3 
e = —p = —nl , 
2 F 2 ' 



(4) 
(5) 



which closes the set of equations for n, v and T. 

We are looking for a set of solutions for a self-similar 
density, flow and temperature profile. Our ansatz for the 



self-similarity means that the non-trivial dependence of 
the solution on the spatial and temporal variables hap- 
pens through the scaling variable 



2 /R 2 (t), 



(6) 



where R(t) is a characteristic, time dependent radial 
scale. Note that the scaling variable x is not to be 
confused with r x , the first component of the space vec- 
tor r = (r x ,r yi r z ). Let us assume rotational symmetry 
and introduce the dimensionless scaling functions v{x), 
(3(x) and T(x) instead of n — n(r,t), v = v(r, t) and 
T = T(r,i), respectively: 



7? 3 

n = f*o-3 v[x), 
R 

v = f v (t)p(x)r, 
T = T f T (t)T(x), 



(7) 

(8) 
(9) 



the time-dependence of the radial scale is written as 

R 2 = R 2 {t) = Rl4>(t), (10) 

and no, Tq and Rq are a constants related to the initial 
conditions. The form of this ansatz is not modified if 
we specify that the minimal "turning point radius" Rq is 
prescribed at a given instant to and the central particle 
density and temperature at this time to is no — n(to,0) 
and To = T(to, 0). This choice of constants implies 



K0)=/3(0) = T(0) = 1, 
/r(*o) = 4>{to) = 1, R a = R{to) 



(11) 
(12) 



The unknown functions fv(t), frit), <p(t) are determined 
as follows. From the continuity equation one obtains that 



fv(t) 



J? 

R ' 



and P(x) 



1. 



(13) 



where ordinary differentiation according to the argument 
is denoted by '. Thus the continuity equation implies 



R ' 



(14) 



regardless of the shape of the density profile. 

Utilizing the equation of state and the continuity equa- 
tion, the energy equation can be rewritten as 



d t T + (vV)T + -TVv 
o 



o. 



(15) 



Inserting the ansatz into this form and making use of 
eq. (Q we obtain the time-dependence of the T as 



1 



Rl 

R 2 



1 



Inserting this to the Euler equation, we obtain 

2T 



f>- 0.5(0') 



mR. 



t + t v - = A 

v 2 



(16) 

(17) 
(18) 



It is trivial to solve these equations for v(x) and 4>{t). 
The solution in terms of physical parameters reads as 



v{x) 



1 



(t~t ) 2 , 



m 



T{x) 



exp 



da 



m(u t y 

2T y r(u) 



(19) 
(20) 



The mean radial flow parameter (ut) controls the asymp- 
totic speed of the expansion, as lim^oo R'(t) = (ut) ■ 
In the above solutions, the functional form of T{x), the 
scaling function of the radial temperature profile, can be 
chosen freely, corresponding to the freedom in the choice 
of the initial boundary conditions. These new hydro so- 
lutions read in terms of physical variables as : 



R 2 {t) = R 2 + (u t ) 2 {t-to) 
v(f,r 



R'{t) 



R(t) 



r = 



(u t ) 2 {t - t ) 



Rl + {u t ) 2 {t-t f V ' 



T(t,r) 
n(t, r) 



T T(x)H = T[t Q ,vR /R(t)} H " 



If 



R\t) 



Rn 



n 



T{t ,o) 



R 3 {t) T[t ,rR /R(t)} 



x exp 



-m(ut) 



2 R 2 {t) 

R 



T{t ,r) 



(21) 
(22) 

(23) 



(24) 



This solution describes a system that contracts for t < to, 
reaches its minimum size Ro at the time of the turning 
point t = to when the flow field vanishes. For later times, 
t > to, the system expands again. The temperature field 
is obtained by rescaling the initial temperature field. 

The above explicit solution for the density profile can 
be written into a simpler form using the scaling variable x 
and the scaling function of the temperature profile, T(x): 



n(t,r) 



no 



Rn 



T(x) R\t) 



exp 



n(u t y 
2T Jo T( 



du 



(25) 



This completes the generic description of our solution of 
NR spherically symmetric fireball hydrodynamics. The 
solution, given by eqs. ( ^l}^5| ), depends on 4 constants 
of integration, Ro, (lit), To and to, as well as on the di- 
mensionless function T(x). Three of the constants of 
integration correspond to the three integrated equations 
(continuity, Euler and energy equation). The constant 



to corresponds to the homogeneity of the equations of 
NR hydrodynamics in time. The free choice of the di- 
mensionless scaling function corresponds to the freedom 
in the specification of the initial boundary conditions. 
Thus we have found infinitely many new solutions of NR 
hydrodynamics, one for each integrable function of one 
variable. Note that each of these solutions is scale invari- 
ant: 



n[t,Y- 



R \ R 3 



RoJ Ro 



n(t ,r), 



T|< ' r !)Ir r(to ' r) ' 



(26) 
(27) 
(28) 



independently of the time t, where R' = R'(to). 

The self-similarity can be utilized to match this solu- 
tion to any given initial time U for a given initial radius 
Ri and a given (rotationally symmetric) temperature dis- 
tribution Tj(r) at time i,. These initial conditions are 



T t (r) =T(U,r), R(t t ) = R t , 



v(r, U) 



R, 



(ui), n(U,0) = n. L . 



This yields the following solutions for <fi and R(t): 



i 



(Ui) 



(t-tj)J + 



{utf 



R{t) 2 = R 2 <t>i(t), T t = Ti(0). 



(29) 
(30) 



(31) 
(32) 



This solution appears instead of eqs. (|19 , 21 ) if the initial 
temperature field, initial radius and initial flow field is 
prescribed instead of the turning point temperature and 
radius. 

The new solution presented above also solves the equa- 
tion of state of an ideal gas for isentropic expansion M : 



[d t + (vV)]pn 



-5/3 _ 



= 0. 



(33) 



Note that the time- dependence of the radius and the so- 
lution for the flow-field is independent of the initial tem- 
perature profile T(to, r). Due to this property of our gen- 
eral solution , we can study various initial conditions in a 
simple manner, assuming various distributions of Tito, r) 
and evaluate its influence on the density profile only. 

Knudsen-gas limit — An essential feature of our so- 
lution is that the central and the asymptotic tempera- 
tures may differ and induce a radial dependence of the 
temperature field. If the initial temperature field is con- 
stant in space, we find a special case of T(x) = 1 and 
T(t a ,r) = T(t ) = T . Eq. © yields that the density 
profile is Gaussian with R 2 G = R 2 To/(m{ut) 2 ): 



N 



[2-KR 2 G {t)] 



3/2 



exp 



2R G {t) 



(34) 



2 



It is straightforward to prove || by evaluating the time 
dependence of the phase-space density, that this case cor- 
responds to a collisionless, free streaming of the particles. 
Hence this parameterization can be applied only to de- 
scribe the post-freeze-out stage of heavy ion collisions. 
Thus, non-trivial hydrodynamical evolution takes place 
in the presence of a non-vanishing temperature gradient 
only, for all rotationally symmetric solutions of the non- 
relativistic hydrodynamical equations that we found. 

The Zimdnyi-Bondorf-Garpman solution — 20 years 
ago, a famous analytic solution of NR hydrodynamics 
was found by Zimanyi, Bondorf and Garpman (ZBG) in 
ref. ||. The temperature and the density for that solu- 
tion explicitly vanishes if |r| is larger than an expanding 
radial scale, |r| > R(t). The ZBG solution has been 
extensively applied to describe particle spectra for NR 
heavy ion collisions. The ZBG solution has the following 
simple initial temperature profile: 



T(t ,r)=T 



1 



R\t) 



if |r| < R{t). 



(35) 



Inserting this to our general solution, eq. (Eq) , we obtain 



n(t, r) = n 



Rq 

R 3 



3 r 



1 



R z 



CV/2-1 



if 



<R(t), (36) 



which corresponds to the ZBG solution with an exponent 
m(u t ) 2 



cy2 - 1 



2T 



1. 



Generalization to d- dimensional expansion - Our ana- 
lytic solution can be generalized to expansions in d di- 
mensions, where d is an arbitrary number. The case 
d = 2 is especially interesting as it corresponds to the 
time evolution of hot and dense hadronic matter in the 
transverse directions in high energy heavy ion ]l|, ^0| , pl| 
and in particle reactions || . When analyzing particle cor- 
relations and spectra in these reactions, it is found that 
the longitudinal expansion is relativistic, (sometimes ex- 
tremely relativistic), however, this longitudinal flow can 
be approximately decoupled from the transversal dynam- 
ics, which is found to be non-relativistic [[TIP]. 

In d dimensions, the divergences of the flow and the 
equation of state are modified as 



Vr d = d, e d = —p. 
The scaling variable is x — r^/R 2 (t), the ansatz is 



(37) 



n = n o TwTtK^) v d = r d f d (t)l3(x). (38) 
R (t) 

The equations for the temperature and the radial scale 
parameter are unchanged. In d dimensions, the energy 
equation eq. ( |l5| ) is modified to 



d t T + (v d V)T + -TVv rf = 0. 
d 



(39) 



From this and the continuity equation we obtain after 
eliminating the singular solutions that 



R'(t) 

R(t) 



(40) 



and the dimension of the expansion, d, cancels from the 
Euler equation. The solutions for d dimensional expan- 
sions are very similar to the d = 3 dimensional expanding 
fireballs, the differences are that in eqs. (^||,^3|) one has 
to utilize d dimensional vectors for coordinate-space and 
flow, and in eq. ( p4j) the density decreases as n cx R(t)~ d . 

The Buda-Lund parameterization — The Buda-Lund 
parameterization (BL) was introduced to describe phe- 
nomenologically the NR transversal flows appearing in 
high energy hadron-hadron and heavy ion reactions, 
given first in refs. |0,^O|. This parameterization keeps 
only the means and the variances of the density dis- 
tribution, the inverse temperature profile and it corre- 
sponds to a scaling transversal flow, such as eq. (p2|). In 
refs. p]j4],|6|,pl| , a radial temperature profile was used, that 
can be rewritten as: 



T(r,r) 



T 



1- 



,AT 



2R l 



1- 



,AT V 



2At 2 



(41) 



where r = (r x ,r y ), the change of the temperature in 
the radial and the temporal directions is controlled by 
the free parameters (AT/T) r and {AT/T) T , respec- 
tively, and C is a constant of normalization. In the 
above parameterization, the longitudinal proper-time, 
t = y/t 2 — rf is introduced instead of the time t , 
corresponding to a relativistic longitudinal expansion. 
The exact Buda-Lund density profile is obtained from 
eqs. (p5pl as 



n(t ,r) = n (1 + (AT/T) r r 2 /(2R 2 Q )) x 



exp 



m(u t y 
2T n 



r / AT. 



Rq 



T 



(42) 



The time-evolution of this density can be obtained from 
simple re-scaling, eq. (p6|), and comparing the BL ansatz 
for the time dependence of the cooling, eq. ([ll|) with the 
exact hydro solution we find 



2Ai 



:(AT/T) r = 



(utf 



(43) 



The exact density distribution in eq. ( |42| ) can be approx- 
imated by a Gaussian, corresponding to the BL Gaus- 
sian ansatz for the density profile in ref. [Q , if terms of 
0(r 4 /R 4 ) can be neglected. (Note that the difference be- 
tween the BL and the ZBG solutions is only C(r 4 /i? 4 ) 
if one sets (AT/T) r = 1 in the BL hydro solution.) The 
BL Gaussian radius is 



1 

R n 



1 

Rn 



m{u t ) 
To 



(AT/T) r > 



(44) 



3 



Hence the BL parameterization is found to correspond 
to a Gaussian approximation to an exact solution of 
the non-relativistic hydrodynamical equations, with con- 
strained values of the cooling parameters (AT/T) T and 
(AT/T) r . If the inequality in eq. (0) is violated, the 



solution for the density profile in eq. (42) will be still 
valid. In this case, the radial density profile will increase 
first before decreasing at large distances. This type of so- 
lution corresponds to a kind of expanding "smoke-ring" 
in the transversal plane, or, to an expanding spherical 
shell if the expansion is three dimensional, as illustrated 
on Figure 1 for R = 5 fm, m{u t } 2 /T = 0.125. For 
(AT/T) r = 0.06 < m{u t ) 2 /T the density has one maxi- 
mum while for (AT/T) r = 0.2 > m(u t ) 2 /T the density 
profile develops two maxima in any transversal direction. 

Application to FOPI data — The FOPI Collaboration 
measured recently the proton-proton correlation func- 
tions at 1.93 AGeV Ni + Ni collisions §]. To interpret 
their data, they utilized, independently, a version of the 
hydrodynamical solution, presented here. They assumed 
a linear flow profile, a Gaussian density distribution and 
a constant temperature. This solution corresponds to a 
collisionless Knudsen gas ||. One expects that a colli- 
sionless approximation breaks down for 1.93 AGeV Ni 
+ Ni collisions. Indeed, only the peak of the FOPI ex- 
perimental proton-proton correlation functions was well 
reproduced by the collisionless model, however, the tails 
had to be excluded from the fit. From the above pre- 
sented solution of the equations of the NR hydrodynam- 
ical systems we know, that the temperature distribution 
must have a radial and a temporal profile if a non-trivial 
hydro evolution is expected. Such a situation was also 
studied in ref. 0], where it was shown that a temper- 
ature profile results in a momentum dependence of the 
effective radius parameters of the proton-proton correla- 
tion functions and in a characteristic bending down of 
the proton- spectra, observable as a deviation from the 
exponential spectrum [^). Hence we predict a deviation 
from the exponential proton spectrum and a momentum 
dependence of the effective source radii, to be observed 
in a more detailed analysis of the FOPI data. Note, that 
comparison with experimental data with the help of a 
hydro solution requires also a condition for the particle 
freeze-out, see ref. || for details. 

Summary — A new family of simple and exact solu- 
tions of the non-relativistic hydrodynamical equations is 
found for rotationally symmetric fireballs, for an ideal 
gas equation of state. Without radial temperature gra- 
dients, the only solution of these equations was found 
to be a collisionless Knudsen gas lacking dynamical evo- 
lution. The results are generalized to expansion in ar- 
bitrary number of spatial dimensions. The Zimanyi- 
Bondorf-Garpman solution Q and the Buda-Lund hy- 
dro parameterization 0] are recovered as special cases. 
For strong enough radial decrease of the temperature, ex- 
panding smoke-ring type of hydro solutions are described. 



Acknowledgments: Exciting discussions are gratefully 
acknowledged with M. Pliimer, J. Balog, P. Csizmadia, 
P. Levai, B. Lukacs, J. Zimanyi and J. Zsigrai. This 
work has been supported by the US - Hungarian Joint 
Fund MAKA 652/1998, by the OTKA grants T026435 
and T024094, the NWO - OTKA grant N 25487 and 
by a former Advanced Research Award of the Fulbright 
Foundation. 



[7 

[§ 

[9 
[10 

[11 



T. Csorgo and B. Lorstad, Phys. Rev. C54 (1996) 1390 
S. Esumi, S. Chapman, H. van Hecke and N. Xu, Phys. 
Rev. C55 (1997) R2163 - 2166 

T. Csorgo, B. L orstad and J. Zim anyi; Phys. Lett. B338 
(1994) 134-140; [nucl-th/9408022; 

J. Helgesson, T. Csorgo, M. Asakawa and B. Lorstad, 
Phys. Rev. C56 (1997) 2626 - 2635 

R. Kotte et al, FOPI Collaboration, Z. Phys. A359 
(1997) 47-54 

N. M. Aga babyan et al, NA22 Collaboration, |hcp 
ex/971 1009| , Phys. Lett. B422 (1998) 359 - 368 



I. G. Bearden et al, NA44 Collaboration, Phys. Rev. Lett. 
78 (1997) 2080-2083 

J. Bondorf, S. Garpman and J. Zimanyi, Nucl. Phys. 
A296 (1978) 320 - 332 



P. Csizmadia, T. Csorgo and B. Lukacs, tmcl-th/9805006 
T. Csorgo and B. Lorstad, Nucl. Phys. A590 (1995) 465c 
T. Csorgo, B. Lorstad, Heavy Ion Phys. 4 (1996) 221-232 



1.4 



1.2 ■ 

1 ■ 

o.s ■ 

0.6 ■ 

0.4 ■ 

0.2 ■ 

- 



<AT7T>,. > m<u,>-/T„ 



<AT/T>,. < m<u.> /T„ 




FIG. 1. Illustration of the development of smoke-ring so- 
lutions for large temperature gradients in exact solutions of 
non-relativistic hydrodynamics. 
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